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Abstract 

We show how nonrelativistic many body techniques can be used to study quantum 
corrections to the classical limit, in particular of the SU{2) Lipkin Model. We show 
that the quantum corrections are essentially of two types: unitary and nonunitary. In 
this work we perform a detailed study of the unitary corrections. They can be cast 
in Hamiltonian form and are shown to double the number of degrees of freedom. As 
a consequence chaotic behavior emerges. We show that this semiquantal chaos is the 
mechanism trough which tunneling is effected. We also show that these corrections 
systematically improve the classical results and propose some quantitative measure of 
this improvement. 
PAC number (21.10) 
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1 Introduction 



One of the most fascinating open problems in the area of dynamical systems is the search 
for an adequate semiclassical description of quantum systems. Several attempts focus on 
obtaining % corrections to the dynamics of the corresponding classical system. A pioneer 
work along this idea is the semiclassical method proposed by Einstein-Brilloin-Keller |jl| 
which leads to quantization rules for integrable systems. Other important contributions are 
approximations to the Feynman path integral formulation used to derive the periodic 
orbit trace formula for chaotic systems: The Gutzwiller trace formula This relates the 
spectrum of quantum systems to a weighted sum over the unstable periodic orbits of the 
classical system. 

Recently semiquantal dynamics has been derived via Ehrenfest's theorem and recast 
as an extended classical system with the fluctuation variables coupled to the average 
variables [Q] j^. A different approach which can be shown to yield identical results is the 
time dependent variational principle where the true solution is approximated by squeezed 
states, the so called gaussian variational approximation ||7|, and Other methods 
for this propose are quoted in ref. |1C]. The Wigner formalism has also recently been 
applied to study diffusive and dissipative type of corrections [11| and [^] . 

The purpose of the present contribution is twofold. The first one is of formal character. 
We show that a mean field expansion in the sense of nonrelativistic many body theories can 
be used as a consistent and systematic tool to analyse the nature of quantum corrections 
to the classical limit. We show that such corrections can be classified in two types: unitary 
quantum corrections, which amount to considering (on top of the dynamics of the average 
values) the dynamics of the width of gaussian wave packets. The number of degrees of 
freedom of this system is therefore doubled, and the semiclassical limit can be cast into the 
form of Hamilton's equations. This unitary dynamics reveals in particular general features 
of quantum kinematics - a centrifugal type barrier involving the width of the packet is the 
classical counterpart of the uncertainty principle - and the resulting dynamics is in general 
chaotic even if we start from a simple integrable system, as we will show. The same result 
in spirit has been obtained in the context of a variational calculation by Pattanyak and 
Schieve Q in the framework of Heisenberg's equations of motion. The second type of 
corrections to the classical limit in the context of many body mean field calculation is non 
unitary in character. It arises from the inclusion of quantum correlation contributions. 
Their inclusion induce, given adequate approximations, a Langevin-type force on top of 
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the Hamiltonian description provided by the unitary evolution. They can thus formally 
account for diffusive and dissipative type behaviors, making immediate contact with the 
work of reference |11]. 

The second purpose of the present work is to explore in detail the unitary type correc- 
tions in the context of the integrable SU{2) Lipkin Model [^]. The reason for choosing 
this particular model, aside its simplicity, resides in the fact that it possesses a well defined 
classical limit when the number of particles goes to infinity [0]. In this case it can 
be rigorously shown that the classical limit corresponds to constructing the Hamiltonian 
dynamics for the parameters of coherent states, a special case of gaussian wave packets. 
In this case we will be safe to perform comparisons between the classical, unitary semi- 
classical and exact dynamics. We show that the derived corrections improve the classical 
results both qualitatively and quantitatively. From the qualitative point of view we show 
that the presence of quantum degrees of freedom (widths of wave packet) induces chaotic 
motion and this is the mechanism through which quantum tunneling is effected, a clear 
quantum correction to the classical description. Also from the qualitative point of view 
we analytically study the time evolution of observable for short times, investigating in 
this way modifications introduced by quantum corrections. From the quantitative point 
of view, in order to assess the time of validity of the approximation we perform a com- 
parative numerical study of the time evolution given by the approximations and the exact 
result for the time evolution of observable. We show that despite of chaotic behavior the 
semiclassical or gaussian approximation gives a better description of the dynamics. We 
also set up a quantitative measure of the time of validity of the approximation accord- 
ing to which the gaussian approximation is better than the classical one and this quality 
increases with N, as it should. 

This paper is divided as follows: section 2 contains the formalism based on which 
we define the semiclassical limit, section 3 contains an application to the SU{2) Lipkin 
Model, section 4 contains a discussion of our numerical results and section 5 contains some 
conclusions. 
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2 The Formalism: 



One of the most widely used method to construct the classical limit of a quantum system 
is by means of coherent states [l^. As is well known such method can be viewed as a 
mean field approximation where the width of the wave packet is minimal and remains 
unchanged during the time evolution. 

The idea of defining the semiclassical limit as some kind of mean field approximation 
is appealing. A very general and unique definition of mean field approximation can be 
given once one requires that the density function which should be used to calculate traces 
be the one which reproduces in exact form all expectation values of one body operators. 
In other words a gaussian density operator jl^] . 

The formalism stated below is a simple application of techniques developed before for 
the treatment of the reduced dynamics of gaussian observable of interacting many boson 



systems in Many Body Nuclear Physics and Quantum Field Theory [17|. We make here 



an option for self-containedness. In the eventual lack of technical details the reader is 



referred to the works in references Ijl^, and [18|. 

Quantum bosonic states are represented by density operators F so that mean values 
of a chosen operator O are given in terms of traces, e.g. 

<0>=Tr{OF) , (1) 
and the basic dynamical equation is the Liouville-von Neumann equation for F 

iF=[H,F] , (2) 

where H is the Hamiltonian of the system. 

Gaussian states are particular densities which are fully determined by the mean values 
of the field operators and their bilinear or quadratic combinations, in addition to statistics. 
For bosonic systems ([a,a^] = 1) the relevant quantities are < > , < a^a > , and 
< a'^a'^ >. The first of these can be conveniently parameterized in terms of two real 
quantities q and p 

<a^>=Tr{a^Fo) = M{q-^^) , (3) 
V 2 

where Fq is the density matrix (it will be defined later) and fiQ is a scale parameter (we set 
liQ = 1). To deal with the other two one can define new Bogoliubov quasiboson operators 
as 
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where fj = rj— < rj >, a = a— < a > and X define the Bogohubov transformation. It is 
given by 

satisfying the normahzation condition 

X^GX = G (6) 

with 

«=(;-i) . «==^ ■ p) 

The preservation of the commutation relations [r],r)^ = 1 requires as usual that the 
transformation coefficients x and y be chosen so that \x\'^ — |yp = 1. This is guaranteed 
by equation (6). 

For the bilinear operators equation (4) implies that 

N = XtRX , (8) 

where 

^ / < jytj) > < 77?7 > \ p_ ^ / < > < aa > \ .g^ 

\^ < fj^fj^ > Krjfj^ > J ' y < a^a} > <aa} > j ' ^ ' 

This Bogoliubov transformation is so chosen that from a direct calculation one has 

TrififiFo) = Tr{f]^fi^Fo) = . (10) 

This gives us also the quantity u =< fj^fj > which is the occupation probability. We can 
thus parameterize < a^a^ > and < a^a > in terms of x, y and u. As a consequence of this 
the dispersion of the quadratures Aq and Ap must depend on x,y and u only 

Aq^ = 7r-[< + >^ -(< a>* + <a >f] = 
2fj,o 

= ^[l+2\y\^-{x*y + y*x)]il + 2,.) , (11) 
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Ap2 = 



-2fio[< {a^ - of > -(< a >* - < a >)^] = 



= 2iso[l + 2\y\^ + {x*y + y*x)]il + 2u) . (12) 

As we will show in the sequel it is possible to define a pair of canonically conjugate 
variables associated to Aq and Ap as follows 

Aq^ = Q^ , 



Q 

The gaussian state so specified is given by 

= i:i">rTT;<Tf^'"<"i ■ 

where ^^T/|n >= >, is the state which contains all the information about the operators 
in question and only this. Notice that the operator ij^fj in the expression of Fq is a 
linear combination of the bosonic operators ,a,a^a,aa,a^a^ and the coefficients are 
functions of q,p,x,y and v. More technical details are given in the didactic work by de 
Cloizeaux [19|. It is a simple matter to check that Fq need not be a pure state and that 



in the limit — > it goes to the pure gaussian state |0 >< 0|. This state precisely 
corresponds to the variational state used in refs. and ^J. The usual coherent state 
approximation corresponds to setting 1^ = and moreover requiring the fullfilment of the 
minimal uncertainty relation, x = 1 and y = for all times. 

Notice that the equation Tr{fjFQ) = Tr{fi^Fo) = is a consequence of the definition of 
fj and therefore imposes no constraints on x and y. Therefore we have four undetermined 
parameters Re(x) ,Im(x), Re(y) and Ini(y) and three conditions to fix them, namely, eq. 
(10) and the normalization condition. The fourth one can be obtained by imposing an 
overall null phase in the state Fq. 

We now sketch the mean field expansion method: In general the state is not of the 
form Fq, but even so can be used to evaluate < >i < ci^o-^ > and < a^a > and hence a 
set of gaussian parameters. In terms of these we can set up a gaussian state Fq and split 
the complete state F as 
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F = Fo + F' (14) 

where Tr(F') = 0, so that F' is a pure correlation part of F. We make essential use of the 
fact that the decomposition (14) can be implemented in terms of a projection operator 
P{t) such that Fo = P{t)F (the explicit form of P{t) for bosons is given in refs. |p and 
[p^ ) to rewrite eq. (14) as a set of coupled equations for Fo and F'{t). This eventually 
allows one to write F'{t) as a function of the past history of the gaussian projection Fo(i): 

F'{t) = F'[Fo{t' < t)] (15) 

Using this result one can close the equations of motion for the gaussian parameters 
again by taking appropriate traces of the Liouville - von Neuman equation (2) and using 
F'{t) = i*b(i) + F'[Po{t' < t)]. The first term will reproduce the mean field result, while 
the second will give rise to additional terms involving memory integrals (refs |Q and pT]]). 
As a matter of fact, the expression for equation (15) is in general not computable even 
in simple model problems without approximations. However a consistent and systematic 
approximation scheme has been constructed based on a criterium of energy conservation 



<H >= Tr{HFo) + Tr{HF') 
so that to every order of approximation one makes sure that d < H > /dt = 0. As 



shown in ref. |2C] and implemented in various systems [|l^] this criterium leads to a 
systematic, controllable expansion around the mean field approximation. 

Since we have defined v = Tr{FQfj^fj) the occupation probability does not evolve on 
time on the mean field approximation level (unitary contributions). We have 

zz> = Tr[ifi), H]Fq = TrH[Fo, fj^j]] , 

and [Fo,fi^ri] = (see eq.(13)) The nonunitary contributions will come from the time 
evolution of quantum correlations, 

z^Tyt^ = ^,>(i) = Tr{[fi'^f], H]F') , (16) 

in the form of an explicit time dependence on the occupation probability ly. This has been 
explicitly studied in the context of the anharmonic oscillator [|l^] and of the Maser Model 
[PU, showing both qualitatively and quantitatively what are the effects of many body 
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correlations at the level of one body observables. In the present work these nonunitary 
contributions have been neglected, since they are of higher order than the unitary ones. It 
is however worthwhile noticing that the inclusion of such corrections may lead to Fokker- 



Planck type equations making thus immediate contact with the works of refs. [11 1 and 



Let us now make the connection between this general mean field expansion and the 
classical limit plus corrections. From the point of view of this formalism the classical limit 
corresponds to the following scheme: 



Exact 



Quantum Mechanics of Many Body Systems 



State Space: 

F=FQ{q,p,Q, P;i>) + F' (correlations) 

Dynamics: 

iF = [H, F] 



Semiclassical 



Most General One Body Mean Field Approximation 



State Space: 

F = Fosc={q,P,Q,P;'y = 0) 



Dynamics: 



Mean Values Dynamics: p = —■iz'Hsc , Q = ^"^s 



Width Dynamics: P = -^'Hsc , Q = ■^p'Hsc 

Hsc = Tr{FoscH) = Hd + Hcorrec 



i 



Classical 



One Body Coherent State Approximation 



Kinematics: 

F = FQ,l{q,p-v = {)) 



Dynamics: 

■ 9 'T_/ ■ 8 



d n_i 



Hci = Tr{FociH) 



where Fq^i is obtained from eq.(13) and the Bogohubov transform (4) by setting x = 1, 
y = and = 0, i.e., FqcI corresponds to a coherent minimum uncertainty state. On 
the other hand -Fqsc incorporates the dynamics of the quadratures < a^a > and < a^a"^ > 
related to the variances (eqs. (11) and (12)) and is a pure state (i^ = 0). The dynamics 
of < a^a > and < a^a^ > enter as a correction {Ticorrec) to the classical one {Tid)- Of 
course considering v ^ Q would enlarge the class of gaussian states so as to encompass 
mixtures. Since there are no rigorous classical results available which cover this generality 
we restrict ourselves to = 0. 

Since on the one body mean field approximation the occupation probability v does 
not depend on time we have an even number of parameters which characterize the time 
evolution of Fq, namely p, x and y - in this case v enter the dynamics as a free parameter. 
This enable one to cast (in a simple way) the variances Aq and Ap in the form of canonical 
variables {Q, P) [jl8| quoted in the above scheme. We define 



X = coshfj + t— 



(17) 



sinh cr 



and then 



(18) 



P 



'(1 + 2Z.) 



(19) 
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In the semiclassical level, defined by -Foso the transformations are the same, setting 
V = 0. The connection of these semiclassical variables with the variables used in refs. Q 
and is made by a simple transformation 

G = Q2 ; n = ^ ;i. = . (21) 

It is a simple matter to check that eqs. (11) and (12), in terms of these variables 
acquire a very simple interpretable form 

Ag2 = Q2 ^ (22) 



Ap^ = P^ + ^q^ . (23) 



It becomes clear how the uncertainty principle will manifest in the semiclassical ap- 
proximation: The fact that the width of the wave packet cannot be zero is translated in 
classical terms by a centrifugal barrier 

^ ■ 

of course this term will be introduced also into the dynamics and coupled to the other 
degrees of freedom. The uncertainty product is given by 



A,Ap=. (1±5£ + P.Q2 > (1±M , (25) 



and the minimal uncertainty situation in the limit ly —>■ is Aq = Ap = We get it 
with P = and Q = . 

The presence of the centrifugal barrier in equation (25) shows that it came from the 
gaussian approximation in a purely kinematical way and therefore does not depend on the 
dynamics H. 

The next step is to obtain the time evolution of the complete set of gaussian parameters. 
This is done by means of the quantum equation of motion 

ij^<0>=Tr{F[0,H]) , (26) 
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where O is any operator we have chosen as relevant. Let us calculate the l.h.s. of eq. (26). 
We calculate the time derivative of the matrix N (eq. (8)) 



— N = — X^RX 28 

dt dt ^ ' 



which gives 



XtRX = N - XtRX - XtRX , (29) 
and with help of the normalization condition X^GX = G, 

XtRX = N - X+GXGN - NGX^GX . (30) 

Writing explicitly the above matrix equation and comparing them element by element, 
we obtain 

ij^fi^'^ =i{xy-xy){l + 2v)=Tr{[fi^fi^,H]Fo) . (31) 

For the condensate < a > the equation of motion is obtained directly from the param- 
eterization (3) and the Bogoliubov transformation as 

at \ z /iQ 

The Hamiltonian dynamics in terms of canonical variables is obtained with transfor- 
mations (17)-(20) which splits eqs. (30) and (31) into a set of four equations for the 
corresponding real quantities q, p, Q and P. 
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3 The Model and the Semiclassical Dynamics: 1/N correc- 
tions 



The SU{2) Lipkin Model ^ 



H = eJ, + y^{Jl + Jl) , (33) 



is characterized in Schwinger's representation by p2| 



H = ^eibh - a^a) + ^ [b^b^aa + a^a^bb] , (34) 
where b^ and are the creation operators for bosons 

[a, = 1 [b, b^] = l . (35) 

Here we have h = 1. The realization of SU(2) algebra in this representation is 

Jz = \ib^b-a^a) , (36) 

J+ = b^a , J_ = a^b. (37) 

The number of particles in the system is given by = b'^b + a'^a and the Casimir 
operator have eigenvalues equal to N/2{N/2 + 1), with N = 2 J. 

In order to obtain the mean field approximation we first define the normalized density 
matrix Fo{a^ a, u^, Vfi), i.e., encompassing mixtures. The corresponding semiclassical 
approximation may be immediately obtained by setting i/j = 0. We have 



where 



and 



Fo = Foa^Fop , (38) 



-Oa = -r^iT^f^ ; Pop = r-^{-r^f^ (39) 



a = Xa{a- < a >) + ya{a'- < >) , (40) 
p = Xb{b- < b >) + Vbib- < b^ >) . (41) 
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The next tedious but straightforward step is to rewrite the Hamiltonian (33) 
in terms of the operators in the equations above {H{a^a,a^a^ ,aa,b^b,b^b\bb) 
H{a''a,a'a\aa,P^P,P^P\PP)) so that the necessary traces are easily calculated. We 
get 

zz>^ = i^< OL^OL >= TrFo[a^a, H] = (42) 

li^O = ^J^< >= TrFoiP'f^, H]=0 (43) 

i^- < OL^OL^ >= ^rFo[a^d^i^] => 
dt ui , J 

^iXaVa - Xaya){^ + 2z^a) = -^XaVai^ + 2l/a)+ (44) 

-Vxlil + 2u^)[xf < /3t >2 +y2 < ^ >2 _xlyb{2up + 2 < /3t X /3 > +1)]+ 
-Vylil + 2ua)[xl < /3 >2 +yf < >2 -x,y;{2u^ + 2 < /jt >< /3 > +1)] , 

lixbVb - Xbyb){l + 2z^/3) = -exbVbil + 2z//3)+ (45) 

-Vxl{l + 2up)[xf < >2 +yl <a>^ -xlya{2ua + 2<a^><a> +!)]+ 

-Vyl{l + 2vp)[xl <a>^ +yf < >^ -Xay*ai2i^a + 2 < >< a > +1)] . 

Prom equations (42) and (43) we see that i^a = i^b = 0- This is directly a consequence 
of the mean field approximation. Thus in the semiclassical level, enters the dynamics 
as a free parameter. For the bosonic condensate the equations of motion read 

ij^<a>=XaTr[a,H]F-yltr[a\H]Fo , (46) 

ij^<b>=XbTr[P,H]F-y;tr[(5\H]Fo , (47) 

^rFo[a^ H] = -]^e[2xaya <a> + y*y„) < >]+ (48) 

-V{yl <a> -yaxl < >)* 
*[xl < /? >2 +yf < /?t >2 -xbyl{2up + 2 < /3t X /3 > +1] 



where 
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-V{xl <a> -XaVl < a"^ >)* 
*[xf < /3t >2 +y2 < ^ >2 -j,*y^^2u(, + 2 < /3t X /5 > +1)] 

TrFo[p\H] = -^€[2xbyb <P> -{xlxb + ylvb) < >]+ (49) 

-V{xl <(5> -XbVl < >)* 
*[xf < >^ +yl < a >^ -xlya{2uc + 2 < >< a> +1)] 

-V{yl <P> -xlyb < >)* 

*[xl <a>^ +yf < >^ -Xay*a{2ua + 2 < >< a> +1)] 

Now the formal definition of the classical limit for type models is in order. 
Prom the SU (2) quasispin operators which characterize the algebra of the model one 
may write ||2^ 

[Ji,Jj\ = —eijkJk (50) 
where Ji = Ji/ J and J = N/2. The classical limit is mathematically defined as 

ymi ^[Ji,Jj\= tijkJk ■ (51) 

rj ^OO % ill 

For finite J the eigenvalues of the ^'s are mapped into the interval [—1,1] and the 
spectra get denser as J increases while h/J decreases. In Schwinger's representation this 



procedure implies that the bosonic operators a and b be scaled as |24] 

"'"-^ ,5.) 

Since we have the Hamiltonian as a function of the Bogoliubov operators 
H{a^a,a^ , aa, (3^ (3, (3(3)) writing down the corresponding equations of motion in 
the canonical variable phase space is straightforward. Firstly we invert the Bogoliubov 
transformation to obtain from parameterization (3) the mean values < a > and < /3 > as 
a function of the scaled quantities = Qi/VJ and = "Pij^fJ {i = a,b). 

< >= \l^[QaiXa + ya) - ip'aiXa - ya)] , (54) 
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<f3^ >=\l^ [q'bixb + Vb) - ip'bixb - Vb)] ■ (55) 

Then we use equations (17)-(20) writing the Bogohubov parameters Xi and yi in terms 
of Qi, Pi and fj. In this way we get 

nsc{qi,Pi,Qi,Pi)=nci{qi,p'd + 
+ jT'^correciQi^Pi^ Pi, = + -j2'^lorrec{Q u Pi\ = ^) , (56) 

Xiliq'a' - P'a')liq'b' - P?) + q'bPWaP'J , 



where 



nLrreM,!/,, Qi, Pi) = ^^{^(Q^ + P^) ' ^iQl + Pi) + - ^^\Qf }+ (^8) 



-X{\{<^b' -p'')[ ^^\q/ - \^Ql - Pi)] + <^MaPa}+ 



-XU{<la-Pa) V ' T^iQt ' P^)] + q'aPaQbPb} , 



niorreciQi, Pi) = ^{[^^ ' '^(QI ' ^D^i^^^^ ' ^(^1 ' + Q^P^QaPa} • 

(59) 

In the above equations x = V{2J) is the scaled interaction parameter. 
In the hmit J — > oo the scaled number of particles 

Af = Tr[{a''a + b'^b)F]/J , 

^f=li<l^'+p^')+li<l^'+p^')+J[^^^^+^^^^^^ , 

(60) 

is a constant of motion. We have 

^^Ar = {n,M} = o , (61) 

with the Poisson brackets defined as 

re/ ^ / ^ r^\-i sr^ df dg dg df Of dg dg df 

{f[qi,Pi,Qi,Pi),9{qi,Pi,Qi,Pi)} = Z^^^ — ^^ + ^7T^~^7T^ • (^^) 

^ dqi dpi dqi dpi dQi dPi dQi dPi 
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Although if J is finite, the Poisson bracket {Tisc^J^} do not vanishes and we have ^ 0. 

As is well known the number of particles N = 2J plays the role of a semiclassical 
expansion parameter in the SU{2) Lipkin Model [^]. Therefore in the limit N ^ oo the 
Hamiltonian T~Lci{qi,Pi) precisely corresponds to the classical limit of the model usually 
taken by means of coherent states . 

The term 'Hlgrredli^Pi^Qii ^i) eq.(56) which is the first order quantum correction 
contains the dynamics of the quadratures coupled to the mean values. Notice that for 
the minimal uncertainty initial condition {g-(0),p^(0), (5i(0) = Pi(0) = 0}, we have 
the total energy equal to the classical one Tisc = T~(-d\t=o {'^lorrec\t=o = 0)- The term 
'^lorreciQii^i) Contains only the dynamics of the quadratures, and if one chose an initial 
condition with the minimal uncertainty situation it does not contribute to the dynamics 
Ti-correcit) — ^ ^he sake of clarity we now reset the variables —>■ qi,Pi)- 

The equations of motion on the {qi,pi,Qi, Pi} phase space are obtained as follows: 
we use (17)-(20) to rewrite the l.h.s. of eqs. (44) and (45) in terms of Qi and Pi and 
parameterization (3) to rewrite the l.h.s. of eqs. (46) and (47) in terms of (ji and pi. 
Doing the same transformations on the r.h.s. of (44)-(47) and comparing their real and 
imaginary parts the two equations for each level a and b split into a set of four real 
equations for the corresponding real quantities qi,Pi, Qi and Pi {i = a,b). 

Qa = -^Pa + Xkai^Y^ + qbPb) " Pa^iql - pI)]+ (63) 



^^^"^ 8JQI 2/^" ^"'^ 



Pa = ^qa - xlPai^Y^ + qbPb) + qa^iqf - Pb)]+ (64) 



8JQl 2J' 



+Xqa[^^-f^-—AQb-Pb)] , 



qb = -^Pb + X[qbi—J— + qaPa) - Pb-iqa - Pa)]+ (65) 

Pb = -^qb - xiPbi^Y^ + qaPa) + qb]^{ql - pI)]+ (66) 
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Qa = -^jPa + - - m+ (67) 



Pa = ^(Qa-^^)+ (68) 

-^(5rpg)[^^ + Qa]-y(^+«.P.)} , 

Qb = +^jPb + - ^(Q^ - )]+ (69) 

+X[y(y^ + 9aPa)-y2(9a-Pa)] ' 

A^^(^^^|^-g.)+ (70) 



2J^ 4Q^ 

2 (U 



1 2 2xr(i+2^^^^i Pb.QaPa^ 
-^(^a -Pa)[^^ + Q.] - y (y^ + 

We next analyze the consequences of these 1/N corrections 
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4 Results: 



What are then the effects of the dynamics of the width (and therefore the uncertainty 
principle) on the classical dynamics of the Lipkin Model? We will discus our results both 
qualitatively and quantitatively. 

Firstly the semiclassical manifestation of the uncertainty principle is the appearance 
of new degrees of freedom whose position Qi cannot be zero. This is formally achieved 
by a repulsive centrifugal-type potential. The inclusion of these new degrees of freedom 
destroys the integrability of the classical dynamics and chaotic behavior emerges. Since 
the Lipkin Model has a regular motion both in the quantum regime and classical limit the 
semiclassical chaos arises as an artifact (a legitimate one) of the approximation. The main 
result we want to stress here is that despite of inducing chaotic behavior quantum (1/iV) 
correction gives a better description of the time evolution of observables. Other interesting 
feature is: Chaos is the mechanism through which quantum properties are effected on the 
semiclassical phase space (such as quantum tunneling effect). 

Let us now describe the classical dynamics in the four dimensional phase space. Beside 
the conservation of energy TCd wc have also the constraint in M. Thus, once it is fixed, 
there exists one, and only one trajectory satisfying J\f =< A'^ > /(J) = 2 with a given 
value of Tici- The existence of these two constants of motion enables one to show in 
the same Poincare section the trajectories for all available energies for a given value of 
the interaction parameter x- In figure 1 we show the well known second order phase 
transition exhibited by the model in its classical limit. For any value of % below the 
critical one |x| < \Xcru\ = 1-0, the invariant tori are all of the same kind and represent 
the rotational aspect of the dynamics (see figure 1(a)). In this case the possible range 
of energies is \Hci\ < 1-0. In figure 1(a) the energies increase from the boundary (where 
Hd = —1.0) to the origin (where Hd = 1-0). For \x\ > \Xcrit\ = 1-0 the possible range 
of energies is enlarged \Hci\ < |(1 + X^)/(2x)|. We still have the rotating trajectories 
\Hci\ < 1-0 (which we label by E^ot in figure 1(b)) and we also have the deformed ones 



-1.0 < Hd < (1 + X^)/(2x) and 1.0 < Hd < -(1 + X^)/(2x) which we label by E, 



and Ejnax respectively in figure 1(b) (x = —6.0). The fixed points of the Poincare map 
associated to the extreme energies are 



= 
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which give us 



2 X 
and 



95 = ; P6 = ±W(1 + -), (71) 



rid = t:" 




g„ = ±y(l + -) ; p„ = 

96 = ; p, = ±J(l--) (72) 




which gives 



11 + x' 



2 X 

The rotational trajectories are isolated from the deformed ones by two separatrices S_ 
(7^cKS-)=-l-0) and S+ {Hci{^+)=1.0). See figure 1(b). 

Introducing the correction terms in Hsc , and therefore the new degrees of freedom 
related to widths {Qi,Pi) coupled to mean values {qi,Pi), the geometrical structure of the 
integrablc system is destroyed and chaotic behavior emerges. The quantity is not a 
constant of motion any more. For increasing J the integrability of the classical limit is 
gradually recovered (see figure 2(a)-(c)) and takes place again only in the limit J oo. 

Of course in the classical domain tunneling effects are completely forbidden. However 
this it is not the case when the quantum corrections we are dealing with are taken into 
account. Let us now look at the low energy orbits (l + x^)/(2x) < "^d < —1-0 (x = —6.0) 
. The classically corresponding invariant tori are localized in the two symmetric regions 
Emm with Pa < oi Pa > exclusively and because of classical integrability the time 
evolution of any chosen initial condition will be confided on its respective region. There- 
fore, destroying integrability is the way the quantum correction works to effect quantum 
tunneling. Choosing an initial condition with energy Hd < —1.0 in the Pa < semiplane 
as an example we show in figure 3(a) its classical Poincare section. For finite values of 
J this same initial condition (evolved semiclassically) is able to access the symmetrical 
region Pa> (figure 3(b)). The quantum observable associated to the transition between 
the semiplanes Pa < and Pa> is = (J+ + J_)/2. Its mean value sign on the Poincare 
section {qt, = , pf, > 0) is determined by the Pa sign. 

Sign[< Jx >] = Sign[^{qbqa + PbPa)] = Sign\pa] . 
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The frequency of such transitions also depends on energy, increasing for energies near 
the separatrix and decreasing as J increase. A quantitative measure of this process may 
be achieved by studying its time scale. We define the confinement time Tc as the time 
interval between two transitions Pa < ^ ^ Pa > ^ ■ We then divide it by the Poincare 
time Tp ,i.e., the amount of time required by starting with an initial condition on the 
Poincare plane and evolving until it reaches the plane again. We take the average over the 
first thousand values of Tc and evaluate the dimensionless quantity Tc/ /Tp as a function of 
J for different energies (see figure 4) . This quantity can also be interpreted as the average 
number of iterations of the Poincare map necessary for a transition to occur. From the 
figure we note that the transitions become scarce as J increase. For any finite value of 
J there must occur a transition, although, for large values of J (or energies close to the 
minimum value) this may require numerically integrating the equations of motion for an 
enormous amount of time. Another interesting feature in figure 4 is the dependence 
on energy. For an energy close to the separatrix (TCd = —1.01) T^/Tp increases slower 
than for lower energies (J~ici = —1.1 and Tid = —1.2). Since tunneling effect manifests 
itself through chaotic motion it must be more conspicuous where chaos (roughly speaking) 
persists, i.e., near the separatrix. 

Notice that here we do not intend to rigorously define a tunneling rate (in terms of 
energy splittings), which is an interesting problem in itself in particular for spin systems, 
we are just characterizing the phenomenon in the gaussian representation. 

An important question that naturally arises in this work concerns the time evolution 
of observables. Since the nonintegrability of the semiclassical description is alien both 
to the quantum and classical dynamics of the system, it is natural to ask whether the 
approximation makes sense quantitatively. We show that despite of introducing chaos the 
inclusion of the quantum degrees of freedom {Qi,Pi) gives a better approximation to the 
time evolution of the observable < Jz > it) which we analyse as an example. 

In figure 5 we display < Jz > {t) (x = —6.0 , J = 4) for the three cases: The exact 
calculation, the semiclassical approximation and the classical result. As can be seen from 
the figure the semiclassical approximation represents an improvement over the classical 
result. We have checked that this is always true for short enough times. The validity 
of the approximation is of course sensitive to the value of x- We next arbitrarily define 
one possible quantitative measure of the accuracy of the approximations. Consider the 
expression 
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approx — pj ~ J ^ TTVnH V'"^j 



J \ Jz ^exact (t) ^ Jz ^approx {t^\dt 
I \ < Jz > exact {t)\dt 

In figure 6 we display for three different values of time, the value of Aapprox as a function 
of 1/J for both the classical as well as for the semiclassical approximation. Notice that the 
error so defined depends linearly on 1/J and the classical calculation lies always above the 
semiclassical one. Figure 6 also shows explicitly that in both cases Aapprox goes to zero as 
J ^ oo. We have also defined a breakdown time in the following way: we fix a maximum 
value for the error A^^^^^, = 0.12 and plot the time when this occurs for several values 
of J in both cases, the classical and the semiclassical (gaussian) approximations. Again, 
according to this measure we see that the gaussian approximation is systematically better, 
i.e., it is valid for longer times (see figure 7). It is interesting to notice that the form of 
the curve is the same for both approximations. The problem of the form of the curve for 



the breakdown time has been the subject of recent investigation |26]. 



21 



5 Conclusions: 



In the present contribution we have shown how a mean field expansion in the sense of 
nonrelativistic many body theories can be used to obtain quantum corrections to the 
classical limit. The unitary time evolution of a gaussian state is shown to contain the clas- 
sical limit plus corrections coming from allowing the width of the wave packet to become 
an independent variable. We have discussed the connection between this approach and 
other approaches in the literature and applied it to the SU{2) Lipkin model. We have 
performed a detailed analysis of the unitary quantum corrections showing that they give 
rise to chaotic behavior, which is essentially the mechanism through which the tunneling 
phenomenon can happen in this context. We have also shown that the quantum correc- 
tions systematically improve the results obtained in the classical limit. The question left 
unexplored in the present work is the effect of nonunitary contributions. We believe this 
is an important next step, i.e., including the time evolution of occupation probabilities in 
the dynamics which is rather natural in the present formalism. It would be interesting to 
cast these contributions in the form of diffusive and dissipative processes. In particular, 
as can be seen from the equations of motion, when t'j 7^ , = the centrifugal barrier 
will be time dependent, affecting thus in particular the tunneling rates. Work along these 
lines is in progress. 

Acknowledgment: We are indebted to Prof. A. F. R. de Toledo Piza for many fruitful 
discussions. 
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Figure Captions: 



Figure 1: 

Poincare section on the plane {qa,Pa) (with = and pb > 0) for the classical 
dynamics 'Hciiqi,Pi) ■{&) for x = -0.5 < Xcrit- (b)for x = -6.0 > Xcrit- See text for 
details. 

Figure 2: 

Poincare section on the plane (qa,Pa) (with = and pb > 0) for the semiclassical 
dynamics Hsciqi, Pi, Qi, Pi) in the case % = —6.0 > Xcrit- Initial conditions for the widths 
are chosen in the minimal uncertainty situation (Qi(0) = -\/l/2, Pi = 0.0). For each initial 
condition the energies are the same as in figure 1(b) Hgc = T^d{t = 0). The values of J 
are: (a) J = 2, (b) J = 8 and (c) J = 12. See text for more details. 

Figure 3: 

Poincare section on the plane {qa,Pa) (with qb = and pb > 0) with x = —6.0 > Xcrit 
for a single initial condition near the separatrix energy TCci{S-) = —1.0 in region 'Emin, 
Pa < 0.0: (a) classical evolution with TCd = —1.1 (in arbitrary e units, e = 1.0). (b) 
semiclassical evolution with Tigc = — 1.1 and J = 9. The initial conditions for the widths 
are set to the minimal uncertainty 

Figure 4: 

Tc/Tp evaluated over the first thousand values of Tc as a function of J for various 
values of Hsc below the classical separatrix energy Hai^-) = —1.0. 

Figure 5: Time evolution of < > (t) {x = —6.0) for the three cases: The exact 
calculation (J = 4), semiclassical approximation (J = 4) and classical result. Time t is 
plotted in arbitrary units. 

Figure 6: 

Error Aapprox evaluated at three different times t (in arbitrary units) and plotted as a 
function of 1/J for both semiclassical and classical approximations, x — —0.5. 

Figure 7: Breakdown time (in arbitrary units) for A^^^^^ = 0.1 plotted as a 
function of 1/J for both semiclassical and classical calculations.x = —0.5. 
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